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ABSTRACT 


In  this  report,  we  present  an  analytical  model  useful  for 
sensor  and  interceptor  trade-off  analysis.  Major  factors  used 
in  this  model  include  sensor  measurement  accuracy,  data  rate, 
interceptor  time  delay  in  responding  to  a  command,  and  inter¬ 
ceptor  control  error  in  executing  a  command.  Guidance  options 
considered  include  command  guidance  and  homing  guidance  whereby 
the  homing  sensor  accuracy  may  either  be  a  constant  or  vary  with 
powers  of  target  range. 
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1. 


INTRODUCTION 


A  guided  interceotor  missile  is  an  essential  component  of  many 
tactical  and  strategic  defense  systems.  Evaluation  of  the  per¬ 
formance  of  an  interceptor  missile  is  usually  a  complicated  matter. 
Before  test  flights  can  be  performed,  exhaustive  computer  simula¬ 
tion  studies  are  usually  conducted  for  use  as  guidelines  on  missile 
designs  and  system  trade-offs.  It  is  often  desired  to  have  a 
simple  analysis  technique  which  can  be  readily  applied  to  provide 
expected  missile  performance  and  preliminary  system  trade-offs 
before  a  sophisticated  and  time  consuming  simulation  program  can 
be  made  available.  It  is  the  purpose  of  this  report  to  discuss 
such  an  analysis  technique. 

The  ultimate  goal  of  an  interceptor  missile  is  to  come  as 
close  to  the  target  vehicle  as  possible  to  insure  that  the  target 
will  be  within  the  lethal  radius  of  the  interceptor  warhead.  The 
closest  approach  between  the  interceptor  and  the  target  is  called  the 
miss  distance.  A  perfect  interceptor  will  ideally  deliver  near 
zero  miss  distances.  An  engagement  system  consists  of  many  mechani¬ 
cal,  electrical,  and  electronic  components;  errors  in  any  components 
will  contribute  to  enlarge  the  miss  distance.  Although  most  of 
these  factors  should  be  included  in  a  detailed  simulation,  they 
cannot  (and  should  not)  be  exactly  represented  in  a  "simple" 
analytical  model.  In  the  model  described  in  tnis  report,  we  in¬ 
clude  the  dominant  system/componer.t  parameters  contributing  to  the 
miss  distance . 
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Dominant  miss  contributors  include  the  guidance  sensor 
measurement  error  and  the  interceptor's  ability  (time  and 
accuracy)  to  respond  to  a  maneuver  command.  A  given  set  of 
parameter  values  characterizing  the  above  variables  will  result 
in  a  miss  distance  value.  On  the  other  hand,  for  a  given  miss 
distance  value,  one  may  find  many  combinations  of  parameter 
values  achieving  the  same  miss  distance.  For  example,  a  poor 
guidance  sensor  can  sometimes  be  compensated  for  with  a  capable 
interceptor  and  vice  versa.  With  this  approach,  one  is  therefore 
able  to  obtain  trade-offs  between  guidance  sensor  and  inter¬ 
ceptor  missile. 

We  emphasize  that  the  work  of  this  report  only  treats  the 
endgame  portion  of  the  engagement.  For  example,  we  assume  that 
the  target  position  uncertainty  (e.g.,  the  handover  basket)  is 
always  contained  in  the  interceptor  reachable  area.  The 
interceptor  is  therefore  not  maneuver  magnitude  limited. 
Discussion  of  interceptor  engagements  from  a  broader  point  of 
view  is  included  in  11, ,  [2] .  The  handover  limit  effect  is 
discussed  in  these  reports. 

This  report  is  organized  as  follows.  In  the  next  section, 
we  give  a  derivation  which  constitutes  the  fundamental  approach 
of  this  analytical  technique.  In  section  3,  we  apply  this 
analysis  to  the  command  guidance  case.  In  section  4,  we  extend 


our  analysis  to  the  homing  guidance  case.  In  a  homing  engagement, 
the  dominant  sensor  error  is  in  the  angle  domain.  The  angle 
measurement  error  can  be  instrumentation  error  limited  (remaining 
constant  as  range  decreases)  or  thermal  noise  limited  (varying 
as  a  power  of  target  range).  All  these  cases  will  be  discussed. 
Finally,  numerical  examples  covering  several  engagement 
scenarios  will  be  given  in  Section  5  to  illustrate  the  method 
and  to  show  how  the  various  error  sources  contribute  to  miss- 
distance. 

An  Appendix,  giving  miss  distance  curves  for  a  variety  of 
systems/component  parameters,  is  included  at  the  end.  This 
appendix  can  be  a  handy  reference  for  system  analysts  evaluating 
interceptor-sensor  trade-off  issues. 
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2 .  FUNDAMENTAL  RESULTS 

In  this  section,  we  present  a  set  of  equations  for  computing  engagement 
miss  distance  given  sensor  measurement  accuracy,  time  interval  between  two 
adjacent  measurements,  missile  time  constant  in  responding  to  a  maneuver 
ccntnand  and  misssile  control  accuracy  in  executing  a  desired  acceleration. 
These  results  can  be  easily  applied  to  ccmmand  guided  interceptors  as 
illustrated  in  the  next  section.  The  methodology  utilized  in  deriving  these 
equations  is  the  same  for  the  honing  guidance  case;  this  will  be  done  in 
the  fourth  section. 

The  interceptor  trajectory  dynamics  are  modeled  as  a  set  of  second 
order  polynomials.  We  will  first  discuss  the  interceptor  tracking  problem 
which  is  then  related  to  miss  distance  calculation.  Limitations  of  this 
analysis  are  also  discussed. 

2.1.  Interceptor  Trajectory  Dynamics 

The  interceptor  trajectory  is  modeled  as  a  second  order  polynomial 
in  each  dimension  with  the  second  derivative  tern  (acceleration)  representing 
the  interceptor  maneuver  control. 

Let  v^  denote  true  interceptor  position  at  time  t^_,  then 

2 


yk  -  P  +  vtk  +  (1/2)  atk 


(2.1) 

where  p,  v,  and  a  are  true  interceptor  position,  velocity  and  acceleration,  at 
t=0,  respectively.  Noise  corrupted  measurements  of  are  represented  by 


zk  =  ^k  +  ?k 

where  is  an  uncorrelated  noise  sequence  with  zero  mean  and  variance 


(2.2) 

2 
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2.2.  The  Interceptor  Position  and  Velocity  Estimates 

The  interceptor  acceleration  a  is  the  command  accelera¬ 
tion  applied  by  the  command  controller.  Suppose  that  the 
acceleration  which  the  controller  wishes  to  apply  to  the  inter¬ 
ceptor  is  an,  the  interceptor  tracking  filter  can  therefore  be 
a  constant  velocity  filter  using  the  following  pseudo-measurements. 


\  = 

zk  ' 

_i_  a  t  2 

2  an  k 

(2.3) 

or 

zk  = 

P  + 

vtk  +  uk 

(2.4) 

uk  = 

+ 

JL(a-a  )  t.  2 

2  n  k 

(2.5) 

where  the  difference,  Aa  =  a- an,  is  the  command  acceleration 

bias  error.  The  pseudo-measurement  noise  sequence  has 

2  2 
mean  1/2  Aatk  and  variance  o z  . 

Assuming  that  the  measurements  are  taken  uniformly  in 

time  and  p,  v,  and  a  correspond  to  the  true  interceptor  dynamics 

at  the  center  of  the  data  interval,  the  following  least  square 

estimates  of  position  and  velocity  can  be  obtained, 


P 

1 

K 

0 

v 

o 

12 

t  2X<k  - 

T  (K-l) K (K+l) — 

where  K  is 


the  total  number  of  pulses  and  T  is  the  time  between 


5 


two  adjacent  measurements.  Using  polynomial  analysis 


we  obtain  ie  following  results: 


(1)  The  position  estimate  £  is  biased  and  the  bias  is 


ilp  =  A  a 


T  ( K — 1 )  (K+l) 
24 


(2.7) 


the  variance  of  p  is 


P 


z 

K 


(2.8) 


(2)  The  velocity  estimate  v  is  unbiased  with  variance 

.  2 


V  T  (K-1)K(K+1) 

A  A  A 

(3)  Using  p,  V,  and  a  ,  the  predicted  position  p  at 
time  t  is  also  biased  and  the  bias  is  p 


(2.9) 


Ap  -  Aa  T  (K-1)(K+1)  +  1/2  Aa  (t  +  ™=iL)2 
P  24  £ 


(2.10) 


The  variance  of  p  is 

P 


2  2 
o_  =  a 

Pp 


i + 


12  ( t  + 


(K-l)T 


2  \ 


T  (K-l) (K+l) 


(2.11) 


where  t=G  corresponds  to  the  terminal  end  of  the  data  interval. 


2.3.  Optimal  Number  of  Interceptor  Tracking  Measurements 

Examining  Eqs.  (2.10)  and  (2.11),  it  is  foun3  that  the  predicted 
position  estimate  bias  increases  with  K  while  its  variance  decreases 
with  K.  An  optimal  K  would  be  the  one  achieving  a  proper  trade-off 


of  these  two  errors.  Assuming  that  the  last  interceptor  command 
which  can  be  effectively  applied  is  limited  by  the  time  between 
measurements  (T)  and  the  interceptor  response  time  constant  (t), 
then  the  optimal  K  is  the  K  minimizing  the  following  terminal 
mean- square  error. 


J  = 


T  (K-l)  (K+l)  . 

—  + 


It  +  T  + 


T  ( K-] 


+  a 


12  (  t  +  T  +  (-—^)T) 

T2 (K-l) (K+l)  i 


(2.12) 


We  note  that  the  interceptor  response  time  (t)  is  the  time  delay 
in  an  interceptor  in  responding  to  a  maneuver  command.  The  partial 
derivative  of  J  with  respect  to  K  is  a  rather  complicated  function 
and  it  is  difficult  to  obtain  a  closed-form  solution  for  K.  Since 
the  total  tracking  time  (K-l)T  is  usually  longer  than  the  inter¬ 
ceptor  response  time  and  the  time  interval  between  measurements , 
we  use  the  following  approximations  in  (2.12)  to  make  the  problem 
tractable 

1)  KT  >>  t 

2)  K  >>  1 


Then  Eq.  (2.12)  becomes 


J  = 


.  2J4  4  o' 

Aa  T  K  +  _ 


36 


(2.12*) 


Taking  the  partial  derivative  of  J  with  respect  to  K  and  sr lving 
for  the  optimal  K  for  minimizing  the  mean-square  error  gives 


where  [  ]  denotes  the  round-off  integer  of  the  enclosed  quantity. 
Notice  that  when  Aa  0,  K  -*•  00 ;  this  implies  that  we  should  use 
all  available  tracking  pulses  when  the  command  bias  is  zero. 

2.4.  Target  Trajectory  Estimation 

The  same  equations  can  be  used  to  determine  the  target 
trajectory  parameters.  In  this  case,  Aa  represents  the  target 
acceleration  error  rather  than  the  interceptor  command  bias.  This 
generally  results  from  imperfect  drag  or  lift  modeling  of  the 
target.  For  the  homing  case,  the  sensor  is  located  on  the 
interceptor;  it  is  therefore  unnecessary  to  separate  tracking 
errors  into  target  end  interceptor  tracks.  The  control  error 
Aa  will  become  the  target  to  interceptor  relative  acceleration 
error  which  encompasses  both  interceptor  control  error  and 
target  drag  and  lift  modeling  errors. 

2.5.  Engagement  Miss  Distance 

In  the  previous  sections,  we  presented  formulas  for 
computing  the  interceptor  tracking  and  prediction  accuracy  and 
determining  the  optimal  number  of  interceptor  tracking  pulses. 

In  this  section,  we  summarize  our  results  for  the  miss  distance 


d 


calculation  as  follows 

1)  Miss  Distance  Bias 


Ax 


Aa 


t  (K-l ) (K+l) 
24 


+  l/2Aa(T  +  T  + 


T  (K-l) 


(2.14) 


(2.14*) 


where  Tj  is  the  interceptor  total  time  interval  defined  in  Eq. 
(2.17)  below 

2)  Miss  Distance  Standard  Deviation 


-  <°/  *  %/>1/2 


(2.15) 


°I  “  az 


1 1  + 


12  (t  +  T  +  -  — ^  -) 
T2(K-1) (K+l) 


1/2 


(2.16) 


*  2c 


1/5 


(2.16*) 


K  = 


6  a . 


2/ 


,  AaT 


(2.17) 


where  and  aTgt  are  interceptor  and  target  tracking  standard 
deviations,  respectively.  As  mentioned  earlier,  the  target 
track  should  have  the  same  characteristics  as  the  interceptor 
track,  we  therefore  do  not  include  an  equation  for  oT  .  Notice 
that  Eqs.  (2.14)  and  (2.16)  are  expressed  in  terms  of  the  raxnber  of 
tracking  pulses  while  Eqs.  (2.14*)  and  (2.16*)  are  in  terms  of 
the  total  tracking  time  interval. 

3)  If  we  use 

x  =  Ax  +  (2.18) 

then  x  is  the  85%  miss  distance,  i.e.,  approximately  85%  of 
the  miss  distances  of  many  random  trials  will  be  less  than  or 
equal  to  x. 

4)  If  we  use 

_____  (2.19) 

x  =  WAx  +  ox 

then  x  is  the  RMS  (Root-Mean-Square)  miss  distance. 

The  above  miss  distance  calculation  formulas  are 
derived  with  the  following  assumptions ; 

a)  The  interceptor  guidance  time  is  sufficiently  long 
and  the  interceptor  launch  is  such  that  the  inter¬ 
ceptor  maneuver  is  capable  of  nulling  the 
predicted  miss. 

b)  The  final  miss  is  the  predicted  position  error 
when  the  time-to-go  is  equal  to  the  sum  of  the 
interceptor  response  time  and  interval  between 
measurements . 
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2.6.  Discussion 

The  above  are  only  some  approximate  formulas  for  quick 
calculation  of  engagement  miss  distances.  It  is  assumed  that  the 
dominant  miss  distance  contributors  are:  1)  position  measure¬ 
ment  error;  2)  interceptor  command  bias;  and  3)  interceptor 
response  time.  Notice  that  these  formulas  are  applicable  to 
command  guided  interceptors  and  homing  interceptors  provided 
that  the  measurement  error  (oz)  is  properly  represented.  When 
applied  to  the  homing  case,  the  position  measurement  error 

(a„)  oecomes  the  relative  error,  and  the  acceleration  error  (Aa) 
z  - 

includes  both  the  command  bias  error  and  the  target  acceleration. 
Shortcomings  of  this  analysis  are: 

1)  It  is  only  a  one-dimensional  analysis.  When  the 
dominant  error  is  indeed  along  a  given  dimension, 
it  may  be  a  close  approximation.  Otherwise,  one 
should  first  calculate  miss  distance  along  each 
independent  axis  then  calculate  the  total  root- 
sum-square  miss  distance  using  these  components 

2)  The  command  acceleration  a  is  calculated  using 
estimated  RV  and  interceptor  states.  This  stochas¬ 
tic  feedback  process  is  not  modeled  explicitly. 

3)  Detailed  dynamic  modeling  of  interceptors  and 
other  effects  are  not  included. 

If  all  the  above  details  are  included,  the  results 
will  not  be  nearly  as  tractable  as  those  presented  above.  Be¬ 
cause  our  model  is  a  simplified  one,  the  resulting  miss  distance 
should  be  treated  as  a  lower  bound  of  actual  performance.  We 
have  compared  the  above  results  with  those  obtained  by  computer 
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simulation  for  several  cases-  The  results  have  always  come  in 
good  agreement. 

Lastly,  we  emphasize  that  the  miss  distance  is  "defined" 
as  the  predicted  target-interceptor  position  error  after  the  last 
effective  maneuver  control.  For  this  reason,  this  model  can 
also  be  used  for  warhead  fuzing  analysis.  In  this  problem,  the 
guidance  sensor  becomes  the  fuze  sensor  and  the  interceptor 
response  time  becomes  the  warhead  reaction  time. 
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3.  APPLICATION  TO  COMMAND  GUIDANCE 


The  equations  derived  in  the  previous  section  are  directly 

applicable  to  command  guidance.  In  this  case,  the  measurement 

error  a  must  be  properly  identified.  Notice  that  in  our  deriva- 
z 

tion  thus  far  we  have  assumed  that  oz  is  a  constant  throughout 
the  engagement. 

A  command  guidance  system  uses  sensor (s)  located 
outside  the  interceptor  vehicle  tracking  both  target  and  inter¬ 
ceptor  and  transmitting  maneuver  commands  to  the  interceptor.  If 
the  tracking  sensor  is  a  radar,  the  dominant  error  component  is 
angle  error.  If  the  target/interceptor  optimal  tracking  interval 
corresponds  to  a  trajectory  length  which  is  much  smaller 
than  the  target/interceptor  range  to  the  tracking  sensor,  then 
the  measurement  error  is  nearly  a  constant  and  can  be  approximated 
with 


a  =  R  0fi 

z  o  fct 


(3.1) 


where  RQ  and  o0  are  mean  range  and  angular  measurement  error 


o 

corresponding  to  the  optimal  tracking  interval. 

If  a  multilateration  tracking  system  is  employed,  the  tracking 
uncertainty  volume  is  the  intersection  of  (at  least)  three 
range  uncertainties.  Let  oR  denote  the  largest  range  uncertainty 
and  3  the  bistatic  angle,  then  the  worst  position  measurement 
uncertainty  is 
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P  <  45 


lVsin  I 

o,  =  !  (3.2) 

(oR/cos  |  p  >  45° 

The  above  formula  should  be  used  as  long  as  the  trilateration 
error  is  smaller  than  the  smallest  cross-range  measurement  error 
made  by  a  monostatic  radar  in  the  multi lateration  tracking  system. 

In  the  numerical  result  section,  we  will  evaluate  miss 
distances  for  a  range  of  oz  values. 
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EXTENSION  TO  HOMING  GUIDANCE 

There  are  four  types  of  homing  sensor  guidance  errors  con¬ 
sidered  in  Ref.  [2],  They  are 


Case 

0 

0-  -  constant 

P 

Glint 

Case 

1 

°p  ”  °eR 

Instrumentation 

Case 

2 

°p  -  Vr2/ra> 

Thermal  (passive  or 
semi-active) 

Case 

3 

Op  -  oa(R3/Ra2> 

Thermal  (active) 

where  o  is  the  position  measurement  uncertainty,  o  is  a  range- 

c  ^ 

invariant  angle  measurement  error  and  is  the  angle  measurement 
error  at  acquisition  (or  a  reference)  range  R^. 

Notice  that  the  case  0  above  is  similar  to  the  command  guid¬ 
ance  case  and  the  equations  of  Section  2  are  readily  applicable, 
we  therefore  will  not  discuss  it  further.  All  remaining  cases 
require  modification  of  equations  of  Section  2. 

The  above  equations  define  an  estimation  problem  with  time- 
varying  measurement  standard  deviations.  One  can  solve  this 
problem  by  applying  the  polynomial  analysis  with  time-varying 
measurement  standard  deviations. 

2 

Let  ok  denote  the  measurement  variance  of  the  k-th  measurement 
?k  (see  (2.3) -(2.5) ) ,  one  can  obtain  the  following  results  for 
the  time-varying  noise  variance  case  after  some  manipulations: 
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(1)  The  position  estimation  bias  at  the  center  of  data 
window  is 


(2)  The  variance  of  the  position  estimate  is 


(4.2) 


(3)  The  variance  of  the  velocity  estimate  is 


(4.3) 


With  the  same  approach  as  in  the  second  section,  the 
optimum  number  of  tracking  pulses  is  the  K  minimizing 

J  -  [Ap  +  ^Aa(T-(^~1)  +T+T)2]2  +  [o  2  +  0  +T+X)2) 


(4.4) 


An  analytical  expression  for  determining  the  optimum  K  is 
difficult  to  obtain.  Since  K  only  takes  integer  values  and  is 
limited  by 


f  1  K  1  2 


(4.5) 


where  is  the  sensor  acquisition  range  and  VR  is  the  relative 
velocity,  an  exhaustive  search  by  stepping  through  its  range  of 
values  appears  to  be  a  reasonable  approach.  The  above  expression 
is  obtained  using  the  fact  that  (1)  the  sensor  needs  at  least 
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two  measurements  to  predict  the  future  trajectory  and  (2)  the 
maximum  number  of  measurements  is  limited  by  the  total  duration 
of  the  engagement . 

The  RMS  miss  distance  is  the  square  root  of  the  minimum  J 
of  (4.4). 

The  evaluation  of  (4.1)-  (4. 3)  will  be  discussed  individually 
in  the  Hollowing  subsections.  In  order  to  simplify  the  manipu¬ 
lations  we  will  use  the  following  approximation  for  replacing  all 
summations : 


^  f  <x,J  *  I  f  (x 

k=l 


k;  *  b-a  J  ^x)dx 


(4.6) 


where  f(-)  denotes  a  function  of  the  enclosed  variable. 

4.1.  instrumentation  Limited  Measurement  Error 
In  this  case,  we  have 

°k  -  Vfl 


(4.7) 


This  equation  implies  that  the  angle  measurement  uncertainty  is 
not  a  function  of  range.  The  target  position  error  is  the  cross 
range  error . 

Examining  (4.1)-(4.3),  one  concludes  that  there  are 
two  specific  summations  to  be  evaluated, 


•i-  Ir 


(4.8] 


k=l  °k 
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(4.9) 


Using  tk  *s  F^/V^  substituting  (4.7)  in  (4.3)  and  (4.9),  and 
applying  (4.6),  one  obtains 


»!  -  -  — 5 — ~ - / 

X  R.  -R  )  J 


% 


dR 


d  <W  1  R  °6  Vb 


S2  *  12  -  ~ 


V<VRa)VR  *R 


/ 


% 


dR  = 


K 


°e  V 


Using  (4.10)  and  (4.11)  in  (4.1)-(4.3)  yields 

2 


AP  = 


AaR  (KT) 

O 

2R7 


2  °6  RaRb 

>  K 

2  °eV 


V  K (KT) 2 


where 


Ra  =  Vr<t+t) 
Rb  =  VR ( t  +KT ) 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


V  =  target-interceptor  relative  velocity. 

K 
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The  optimum  K  is  the  K  which  minimizes  the  J  of  (4.4) 
with  terms  defined  above;  the  miss  distance  is  the  square  root  of 


4.2.  Thermal  Error  Limited  Case  -  Passive  and 
Semi-Active  Homing 

In  this  case,  the  angle  measurement  standard  deviation 
is  inversely  proportional  to  the  range  to  the  target  through 
the  relationship  with  signal-to-noise  ratio.  Let  denote  the 
angle  measurement  accuracy  at  a  reference  range  ,  one  then  has 
the  cross  range  measurement  error  as 

2 

°k  “  °A  ST  (4'17) 

f 

Applying  the  same  derivation,  we  obtain  these  following 


(4.18) 


dR  K  RA2(RbH3) 

n 


R’  °A2 (Rb_Ra) 3Ra3Rb3 


K  R. 


c‘(R_  — R  )V‘  J0  R2  o  2r  r  vd2 
A  d  a  R  R_  AabR 

a 


(4.19) 


3lRb-Ra)Ra3V 


(4.20) 


(4.21) 


°A  Vb 
- _  - 

K  Ra  (KT) 


(4.22) 


Similarly,  one  applies  (4 . 20) - (4 . 22)  in  (4.4)  and 
searches  for  the  minimum  to  obtain  the  miss  distance. 


4.3.  Thermal  Error  Limited  Case  -  Active  Homing 

In  active  homing,  the  square-root  of  the  signal-to- 
noise  ratio  is  inversely  proportional  to  the  range  squared.  The 
cross  range  error  is  therefore 

,  3 


°k  =  °A  fl 

ra 


(4.23) 


Similarly,  one  obtains  the  following  results: 


I,  = 


<Rb 


V  A  dR. 

b-v  i  r6 


K  Raj<Rb5-Ra^ 


= 


(Rb 


Ap  - 


rr4  A 

'VVR2  i 


(4.24) 


dR.  K  Ra4(Rb3-Ra3> 
r4  “a2 (Rb'Ra’ 3Ra3Rb3vR2 


ia(KT) 


5(Rb3-Ra3)Ra2 


(4.25) 

(4.26) 


2  _  a2 

K  R  4 
A 


(4.27) 


20 


r 


K  R*  (KT) 
A 


5 


3(VRa)Ra  Rb 

Rb3RaT 


(4.28) 


The  miss  distance  is  the  square  root  of  the.  minimum  J 
of  (4.4)  with  terms  defined  above. 


4.4.  Summary 

in  this  section,  we  have  extended  the  miss  distance 
equations  of  Section  2  to  the  homing  guidance  case.  Equations  for 
evaluating  estimation  bias  and  variances  are  given.  An  integral 
for  approximating  summations  is  used  in  deriving  these  results. 

The  root-mean-square  miss  distance  is  obtained  with  an  exhaustive 
search  of  the  optimum  number  of  measurements  minimizing  the 
terminal  miss  distance. 

We  re-emphasize  that  when  these  equations  are  applied 
for  the  homing  case,  Aa  is  the  target  to  interceptor  relative 
acceleration  uncertainty  and  it  will  no  longer  be  necessary  to 
separate  miss  distance  into  target  track  and  interceptor  track 
individually. 


5. 


NUMERICAL  RESULTS 


In  this  section,  we  present  numerical  results  illustrating 
system/component  parameter  trade-offs.  The  results  are  in  terms 
of  miss  distance  vs.  measurement  error  (o  or  oy)  with  interceptor 
response  time  t  as  a  parameter. 

In  Figs.  5. 1-5. 4,  we  present  several  comparisons  for  various 

system  and  design  parameters.  A  comparison  of  various  guidance 

modes  is  given  in  Fig.  5.1.  These  are  for  RA=10km,  VR=10km/s, 

2 

control  acceleration  error  of  1  m/s  ,  and  pulse  repetition  rate  of 
50  (T= . 02  s) .  All  pertinent  parameter  values  are  shown  on  the 
upper  right-hand  corner  of  the  figures.  A  set  of  nominal  values 
used  for  the  interceptor  response  time  is  shown  on  the  upper 
left-hand  corner  of  the  figures.  Notice  that  there  is  an  upper 
limit  on  the  interceptor  response  time  (Tmax) .  It  is  limited  by 
the  total  engagment  time  (=R  /V„)  and  the  requirement  of 

A  K 

processing  at  least  two  measurements  for  trajectory  prediction. 

The  bottom  curve  corresponds  to  the  smallest  response  time 
indicated  in  the  upper  left-hand  corner  and  the  upper  curve  corres¬ 
ponds  to  the  largest  response  time  used  for  this  case.  The 
response  time  for  both  curves  are  labelled.  Response  times  for 
intermediate  curves  follow  those  shown  in  the  upper  left-hand 
corner.  One  can  compare  parameter  values  for  achieving  a  giving 
miss  distance  using  these  results.  Notice  that  for  the  command 
guidance  case,  a  faster  interceptor  cannot  compensate  for  poorer 
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measurement  accuracy  when  i  is  less  than  about  .3  second.  On  the 
contrary,  a  responsive  interceptor  can  compensate  for  a  poor 
sensor  for  the  homing  case.  Especially  for  the  thermal  noise 
limited  cases  very  poor  sensor  accuracy  can  be  tolerated  with  a 
good  interceptor.  This  is  because  the  cross  range  measurement 
error  is  proportional  to  target  range  and  sensor  angle  measurement 
error.  For  the  thermal  noise  limited  cases,  the  angle  measurement 
error  improves  significantly  with  decreasing  range.  When  target 
and  interceptor  are  ery  close,  the  cross  range  measurement  error 
becomes  very  small  even  for  large  angle  measurement  error  at 
acquisition. 

One  must  keep  in  mind  that  any  practical  homing  sensor  will 
have  glint,  instrumentation  and  thermal  errors  which  nu  st  be 
combined  to  obtain  the  overall  sensor  accuracy. 

Comparing  Fig.  5 . 1 (b) -5 . 1 (d) ,  it  is  seen  that  the  instrumen¬ 
tation  limited  homing  case  gives  the  most  strict  requirement  on 
both  sensor  and  interceptor.  We  will  therefore  concentrate  on 
the  instrumentation  limited  homing  case  for  the  remainder  of  the 
report. 

In  Fig.  5.2,  we  examine  the  effect  of  closing  velocity  for 
2 

ra=10  km,  &a=l  m/s  and  T=.02s.  The  miss  distance  decreases  with 
decreasing  closing  velocity.  It  is  not  a  strong  function  of  t 
until  t  becomes  comparable  with  the  time-to-go  (RA/VR) •  Miss 
distance  is  reduced  almost  linearly  with  sensor  angular  error 
except  for  very  small  errors. 


23 


towawsiSSSH 


STANCE 


INSTRUMENTATION  LIMITED 

r  01  ,  .02,  .05.  .  1  ,  .2,  .5, 
1..2..5..10. 


-  mr  ad 


Fig.  5.2(c)  Effect  of  initial  time-to-go 


The  effect  of  data  rate  is  shown  in  Fig.  5.3.  These  results 

are  for  RA=10  km,  Vc=10  *tIn/,s'  Aa=0.1  m/s  ,  and  data  rates  of 

20,  50,  and  100  measurements  per  second.  There  is  a  uniform 

reduction  of  miss  distance  at  higher  data  rates.  Figure  5.4 

gives  the  dependence  on  control  error.  For  large  o  the  miss- 

distance  is  independent  of  Aa  but  as  crQ  decreases,  the  dependence 

on  Aa  is  almost  linear.  These  results  are  for  RA=10  km,  Vc=3  km/s, 

2 

T=0.01  second,  and  control  error  of  .1,  1  ,  and  10  m/s  . 

A  cross-plot  of  miss-distance  as  a  function  of  missile  response 

time,  t,  with  as  a  parameter  is  shown  in  Fig.  5.5.  Note 

that  the  dependence  is  approximately  /T  for  small  t  and  large 

oM .  in  addition,  it  shows  the  effect  of  t  which  is  determined 
0  max 

by  the  available  engagement  time. 
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6. 


SUMMARY 


In  this  report,  we  have  presented  an  analytical  model  use¬ 
ful  for  sensor  and  interceptor  trade-off  analysis.  This  model 
can  treat  various  guidance  modes  such  as  command  guidance  and 
homing  guidance  whereby  the  homing  sensor  accuracy  may  either 
be  a  constant  or  vary  with  powers  of  target  range. 

The  advantage  of  this  model  is  its  simplicity.  Although 
simplified,  it  does  include  major  miss  contributors  including 
sensor  accuracy,  data  rate,  interceptor  time  delay  in  responding 
to  a  given  command,  and  interceptor  control  execution  error. 

The  results  of  this  model  therefore  give  a  tight  lower  bound 
to  the  actual  performance. 

This  model  also  assumes  that  the  interceptor  is  not  maneuver 
force  limited  and  has  sufficient  fuel  for  nulling  the  initial 
miss.  This  requirement  corresponds  to  the  fact  that  the  pre¬ 
guidance  sensor  provides  good  handover,  thereby  eliminating  the 
need  of  excessively  high  acceleration  maneuvers. 

This  model  can  also  be  used  for  warhead  fuzing  analysis. 

In  this  application,  one  replaces  the  guidance  sensor  with  the 
fuze  sensor  and  the  interceptor  response  time  with  the  warhead 
response  time. 

An  appendix  containing  miss  distance  results  for  a  wide 
range  of  parameters  is  also  given. 
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We  summarize  some  approximate  scaling  rules  for  miss-distance 


as  a  function  of  various 
Sensor  Error: 

Missile  Response 
Time : 

Data  Rate: 

Closing  Velocity: 
Command  Bias: 


engagement  parameters  below: 

Miss  distance  ~  unless  is  very 
small 

Miss  distance  ~  /t"  unless  t  is 
comparable  to  R  /V0 

Miss  distance  ~  /t  unless  0^  is  very 
small 

Miss  distance  —  V_  unless  i  is 
comparable  to  RA/vR 

Miss  distance  —  Aa  when  is  very 
small 
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APPENDIX  A 


A  Collejtion  of  Miss  Distance  Results  for 
Homing  Guidance  with  Instrumentation  Error  Limited  Sensor 

In  this  appendix,  we  give  a  set  oi  miss  distance  curves 
covering  a  wide  range  of  design  parameters  for  the  homing  guidance 
instrumentation  error  limited  case.  An  index  to  the  figure 
numbers  and  parameters  used  is  given  in  Table  A.l. 

We  re-emphasize  that  these  results  apply  for  the  case  when 
the  interceptor  is  not  maneuver  force  limited. 


4b 


INDEX  ON  CASES  EXAMINED 


M 

m 

00 

H  1 

•H 

CN 

o  1 

• 

* 

* 

■ 

■ 

< 

< 

< 

r- 

VO 

H 

CO 

CN 

O 

• 

' 

< 

< 

< 

H 

VO 

tn 

m 

*— i 

CN 

o 

• 

a 

a 

• 

< 

< 

< 

o 

rH 

in 

*h 

i c 

• — i 

CN 

o 

• 

• 

» 

• 

< 

< 

< 

in 

m 

m 

o 

• 

•h 

cn 

* 

< 

• 

< 

< 

m 

cn 

i—i 

• 

rH 

<N 

o 

< 

a 

< 

< 

CN 

rH 

in 

m 

*— 1 

CN 

o 

# 

• 

< 

c 

in 

O 

» — i 

CN 

r— ( 

CN 

o 

• 

. 

. 

< 

< 

< 

o 

CXi 

in 

•h 

«H 

o 

• 

• 

< 

< 

< 

Ifl 

1*0- 

0  *  T 

0  *01 

£ 

u 

£ 

0 

1 

cn 

X 

> 

1 

Eh 

2S/“ 

-  Joaaa  pueuiuioa 

46 


Note:  (1)  Entries  are  figure  numbers. 

(2)  All  cases  are  homing  with  instrumentation  limit. 


.fJUt . 


ACKNOWLEDGMENTS 


The  author  would  like  to  express  his  thanks  to  Drs.  R.  W. 
Miller  and  S.  D.  Weiner  who  have  provided  motivating  and  stimu¬ 
lating  discussions  for  this  work.  Dr.  Weiner  has  also  provided 
careful  review  of  the  manuscript  with  suggestions  which  have 
greatly  improved  the  quality  of  this  report.  The  skillful  typing 
and  preparation  of  the  manuscript  by  Lynne  Taibbi  is  also 
appreciated . 


REFERENCES 


( 1)  S.  D.  Weiner,  "A  Simple  Graphical  Model  for  Analyzing  Radar 
Homing  Interceptor  Engagements" , Technical  Note  1979~82  , 
Lincoln  Laboratory,  M.I.T.  (17  December  1979),  AD-A082460/ 7 . 

(2)  S.  D.  Weiner,  "A  Model  of  Interceptor  Engagements  Including 
Warhead  and  Fuze  Interactions  ", Technical  Note  J.  98  0-4  5  > 
Lincoln  Laboratory,  M.I.T.  (18  September  1980),  AD-A093174/1 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (then  Dot a  Entered) 


REPORT  DOCUMENTATION  PACE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  REPORT  NUMBER 


i.  GOVT  ACCESSION  NO.  I  3.  RECIPIENT'S  CATALOG  NUMBER 


ESD-TR-81-365 


4-  TITLE  (and  Subtitle) 


AH. 


3  C  Wfc 


A  Model  for  Sensor-Interceptor  Trade-Off  Analysis 


7.  AUTHOR!!/ 


Chaw-BIng  Chang 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 
Lincoln  Laboratory,  M.I.T. 

P.O.  Box  73 
Lexington,  MA  02173 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Ballistic  Missile  Defense  Program  Office 
Department  of  the  Army 
5001  Eisenhower  Avenue 
Alexantl  .a,  VA  22333 

14.  MONITORING  AGENCY  NAME  A  ADDRESS  Fi'  different  from  Controlling  Office/ 

Electronic  Systems  Division 
Hanscom  AFB,  MA  01731 


14.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  enured  in  Black  20,  if  different  from  Report) 


II.  SUPPLEMENTARY  NOYES 


tl.  KEY  WORDS  F Continue  on  reverse  tide  if  necessary  and  identify  by  black  number J 

sensor-interceptor  trade-off  analysis 
data  rate 

interceptor  control  error 


5.  TYPE  OF  REPORT  I  PERIOD  COVERED 

Technical  Report 

k.  PERFORMING  ORG.  REPORT  NUMBER 

Technical  Report  599 

l.  CONTRACT  OR  GRANT  NUMBERFi/ 


F19628-80-C-0002 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  l  WORK  UNIT  NUMBERS 

Program  Element  Nos.63304A 
and  63308A 

12.  REPORT  DATE 

18  January  1982 

13.  number  of  pages 

82 _ 

IS.  SECURITY  CLASS,  (of  tbit  report) 

Unclassified 


15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


sensor  measurement  accuracy 
interceptor  time  delay 
homing  guidance 


70.  ABSTRACT  (Continue  on  reverse  tide  if  ne  cess  my  nod  identify  by  block  n'unber) 


In  this  report,  we  present  an  analytical  model  useful  for  sensor  and  interceptor  trade-off  analysis. 
Major  factors  used  ir.  this  model  include  sensor  measurement  accuracy,  data  rate,  interceptor  time 
delay  in  responding  to  a  command,  and  interceptor  control  error  in  executing  a  command.  Guidance 
options  considered  include  command  guidance  and  homing  guidance  whereby  the  homing  sensor  accuracy 
may  either  be  a  constant  or  vary  with  powers  of  target  range. 


DD  F0RM  U73 

1  JAN  73 


EDITION  OF  I  NOV  65  IS  OBSOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  FACE  (When  Data  Enter*  Jj 


-"-rr 


